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ABSTRACT 

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations 
by relying on the emission-weighted temperature Tew, instead of on the spectroscopic X-ray temperature Tspec, 
which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is 
thermally complex, Tew fails at reproducing Tspec, and propose a new formula, the spectroscopic-like temperature, 
Tsi, which approximates Tspec better than a few per cent. By analyzing a set of hydrodynamical simulations 
of galaxy clusters, we find that is lower than Tgw by 20-30 per cent. As a consequence, the normalization 
of the M-Tsi relation from the simulations is larger than the observed one by about 50 per cent. If masses in 
simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and /3-model gas 
density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and 
significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer 
the amplitude of the power spectrum from the X-ray temperature function could bias low erg by 10-20 per cent. 
This may alleviate the tension between the value of dg inferred from the cluster number density and those from 
cosmic microwave background and large scale structure. 

Subject headings: galaxies: clusters: general - X-rays: galaxies masses - cosmology: observations 



1. INTRODUCTION 

Thanks to their position at the top of the hierarchy of cosmic 
structure formation, galaxy clusters represent powerful tools to 
constrain cosmological parameters, like the present value of 
the matter density parameter fiom, and the normalization of the 
primordial power spectrum, usually expressed in terms of erg, 
the r.m.s. matter fluctuation in a sphere of radius of 8/i~' Mpc 
(Rosati et al. 2002; Voit 2004 for recent reviews). The funda- 
mental quantity entering theoretical models is the cluster mass, 
but from the observational point of view its high-precision mea- 
surement is still precluded by the presence of different system- 
atic effects affecting the lensing data and the dynamical models. 
A possible short cut to overcome this difficulty is represented 
by the determination of the cluster temperature from X-ray ob- 
servations. In fact simple theoretical arguments, supported by 
numerical hydro-N-body simulations, suggest the existence for 
virialized gravitational systems of a tight relation between mass 
M and temperature T. 

The enormous progresses in the angular and spectral reso- 
lutions of the latest generation X-ray satellites, Chandra and 
XMM-Newton, allowed us to obtain temperature measurements 
for extended sets of clusters, showing in many cases the pres- 
ence of a rather complex thermal structure. The theoretical 
interpretation of these high-quality data requires state-of-the 
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art numerical simulations, in which the physical processes act- 
ing on the gas component (radiative cooling, non-gravitational 
heating, etc.) are faithfully reproduced. The comparison be- 
tween observations and simulations is, however, not always 
simple. As for X-ray data, the cluster gas temperature is ob- 
tained by fitting the source spectrum to a suitable thermal model. 
This procedure cannot be reproduced in the analysis of simula- 
tions, without using devoted software packages (like X-MAS; 
Gardini et al. 2004), which simulate the spectral properties of 
the simulated clusters taking into account the instrument re- 
sponse, as well as sky and instrumental backgrounds. To over- 
come this problem, it has been of common use to resort to the 
emission-weighted temperature, 

_ jAjTWTdV 
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Here n is the gas density and A(r) is the cooling function, for 
which most of the works in the literature use A(r) cx VT, as- 
suming bremsstrahlung emission. 

Recently, Mazzotta et al. (2004, M04 hereafter) demonstrated 
that, if clusters have a complex thermal structure, Tgw always 
overestimates the spectroscopic temperature Tspec obtained from 

spec 



X-ray observations. The discrepancy between Tgw and 



strongly depends on the complexity of the cluster thermal struc 
ture (see also Mathiesen & Evrard 2001; Gardini et al. 2004). 
M04 proposed a new definition of temperature to be imple- 
mented in the analysis of simulations, the spectroscopic-like 
temperature, 

" Jn^T'^/TV^dV ' 

When applied to simulated spectra of clusters hotter than 2-3 
keV observed by Chandra or XMM-Newton , this equation, 
with a = 0.75, yields a temperature within a few per cent from 
Tspec- M04 also showed that this result is almost independent of 
the details of the detector response function, as long as its spec- 
tral properties are similar to or worse than the Chandra ones. 
In this way, Eq. 2 can be reliably applied to the detectors on 
board of Beppo-SAX and ASCA. 
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In this Letter we address the impHcations of using Ts\ on the 
mass-temperature relation and on the X-ray temperature func- 
tion (XTF) for an extended set of simulated clusters. These 
results are then used to discuss the effect on the estimate of 
from the XTF. 

2. THE SIMULATED CLUSTERS 

The sample of simulated galaxy clusters used in this paper 
has been extracted by the large-scale cosmological hydro-N- 
body simulation of a "concordance" ACDM model (fiom = 0.3, 
riQA = 0.7, f^ob = 0.019/!"^ h = 0.7, (78 = 0.8), presented in Bor- 
gani et al. (2004, B04 hereafter). Here we give only a short 
summary of its characteristics, and we refer to that paper for 
more details. The run, made by using the massively parallel 
Tree-nSPH code P-GADGET2 (Springel et al. 2001; Springel 
2004, in preparation), follows the evolution of 480^ dark mat- 
ter particles and an equal number of gas particles in a peri- 
odic cube of size \92h~^ Mpc. The mass of the gas particles 
is nigas = 6.89 x 1O^/!~'M0, while the the Plummer-equivalent 
force softening is 7.5/1"' kpc at z = 0. Besides gravity and hy- 
drodynamics, the simulation includes the treatment of radia- 
tive cooling, the effect of a uniform time-dependent UV back- 
ground, a sub-resolution model for star formation from a mul- 
tiphase interstellar medium, as well as galactic winds powered 
by SN explosions (Springel & Hernquist 2003). At z = we ex- 
tract a set of 95 temperature-selected clusters (hereafter Set#l) 
with Tew > 2 keV, value below which the single spectroscopic 
temperature for thermally complex systems cannot be defined 
(see discussion in M04). 

Due to the finite box-size, the largest cluster found in the 
cosmological simulation has Jew ~ 7 keV. In order to extend our 
analysis to more massive and hotter systems, as required for the 
study of the M-T relation, we further include 4 more galaxy 
clusters (hereafter Set#2) having Mvir > IQ^^H'^Mq and be- 
longing to a different set of hydro-N-body simulations (Dolag 
et al. 2004 in preparation). Since these objects have been 
obtained by re-simulating at high resolution a patch of a pre- 
existing cosmological simulation, they have a better mass res- 
olution, with mgas = 1.69 x 10^/!~'Mq. These simulations have 
been performed by using the same code with the same choice 
of the parameters defining star-formation and feedback. 

Therefore, our total sample (Set#l and #2) comprises 99 ob- 
jects spanning the temperature range 2 "x.7e„ 'x.13 keV. 

3. RESULTS 

For each cluster in our sample we calculate both the emission- 
weighted and the spectroscopic-like temperature. As for Tew, it 
is computed as in B04 by applying Eq. 1, with the cooling func- 
tion estimated in the [0.5-10] keV energy band so as to account 
for the typical energy coverage of X-ray detectors. For simi- 
lar reasons, we apply Eq. 2 to compute T^x by removing the gas 
particles with T < 0.5 keV. We note that, since the mean tem- 
perature of each cluster here considered is higher than 2 keV, 
the excluded particles usually provide a small fraction of the 
total emission measure. Therefore, this procedure does not sig- 
nificantly affect the values obtained for both Tew and T^i. 

In Fig. 1 we plot the value of T^i versus Tew, which confirms 
the result of M04, that for all clusters Tsi is smaller than T^^. 
The distribution of the points in Fig. 1 suggests that the two es- 
timators can be roughly related by a linear relation. By applying 
a fit, we find the relation, 

T,i = (0.70 ± 0.01)rew+ (0.29 ± 0.05) , (3) 
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Fig. 1. — The comparison between emission-weighted temperature, Tew, 
and spectroscopic-like temperature, Ts\. Circles and triangles refer to objects 
from Set#l and Set#2, respectively. The dotted line is the best-fit linear re- 
gression (see Eq. 3), while the dashed line corresponds to r,i = Tew 



shown as dotted line in Fig. 1. Therefore, for objects with tem- 
peratures of about 6 keV, the disagreement can be as large as 2 
keV. As also discussed by M04, the relative mismatch between 
Tew and T^j depends only on the thermal structure of the ICM 
and, in particular, increases as a function of its complexity. For 
this reason, we expect Eq. 3 to hold for the physics considered 
in our simulation sets. Any physical process (e.g. stronger feed- 
back or thermal conduction), which produces a more isothermal 
ICM, is expected to provide a smaller difference between the 
two definitions of temperature. 

The existence of a systematic trend between the two differ- 
ent estimates of the cluster temperature has interesting conse- 
quences on the scaling relation between mass and temperature. 
Adiabatic hydrodynamical simulations (e.g. Evrard et al. 1996; 
Bryan & Norman 1998) find a normalization which is about 
40 per cent higher than observed (e.g. Horner et al. 1999; 
Nevalainen et al. 2001; Finoguenov et al. 2001; Allen et al. 
2001). This result has suggested the necessity of including in 
numerical simulations the treatment of extra physical processes 
for the gas component, like non-gravitational heating and ra- 
diative cooling. For instance, B04 showed that considering both 
effects can be efficient in reducing the normalization of the sim- 
ulated M-T relation, which however turned out to be still 20 
per cent higher than observed. In that paper, as in most nu- 
merical works on the X-ray properties of galaxy clusters, the 
emission-weighted temperature is used to construct the M-T 
relation. However, since Tew tends to overestimate Tsi by a fair 
amount, we expect the actual M-Tsi relation from simulations to 
be even more discrepant with respect to observations. This dis- 
crepancy is indeed apparent in the left panel of Figure 2, where 
the simulated M-T relation at overdensity p/pcr = 500' is com- 
pared to the observational results by Finoguenov et al. (2001), 
who measured temperatures from ASCA data. By modelling 
the scaling relation with 
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' Here per is the cosmological critical density. 



M5oo=Mo 



(4) 



3 




Fig. 2. — The mass-temperature relation at p/pcr = 500, in simulations 
(filled circles and triangles for Set#l and Set#2, respectively) and for the ob- 
servational data by Finoguenov et al. (2001, squares with errorbars). The left 
panel is for the true masses of simulated clusters; the right panel is for masses 
of simulated clusters estimated by adopting the same procedure applied by 
Finoguenov et al. to observational data (see text). 



a log-log least-square fitting for the simulated clusters gives 
a = 1.66 ±0.09 and log(Mo//i-'MQ) = 13.54 ±0.03, while 
log(Mo//!"'M0) = 13.60±0.07by fixing the slope to the value 
a = 1 .5 expected from the condition of hydrostatic equilibrium. 

While masses of simulated clusters have been computed here 
by summing over all the particles, masses of X-ray observed 
clusters are usually estimated from data (e.g., Markevitch 1998; 
Nevalainen et al. 2000; Finoguenov et al. 2001) by applying 
the condition of hydrostatic equilibrium to a spherical gas dis- 
tribution described by a /3-model (Cavaliere & Fusco-Femiano 
1976), with the equation of state having the polytropic form. In 
this way, the total self-gravitating mass within the radius r is 

M«r) = 1.11 X 10^Vfit7 ?'^'? .f , r'Mo. (5) 

keV h 'Mpc l+x^ 

Here T{r) is the temperature at r, /3fit is the fitted slope of the 
gas density profile, x = r/rc is the radial coordinate in units of 
the core radius r^, and 7 is the effective polytropic index. 

A number of analyses of hydrodynamical simulations of clus- 
ters (e.g., Muanwong et al. 2002; B04; Rasia et al. 2004; Kay 
et al. 2004; cf. also Ascasibar et al. 2003) have shown that 
Eq.5 underestimates the actual cluster mass by about 20 per 
cent. Here we repeat the same analysis as in B04, but using 
Tsi in Eq.5 instead of Tgw. Although the resulting M-T rela- 
tion (shown in the right panel of Figure 2) is still slightly higher 
than the observed one, we find that the amplitude of the M-T 
relation decreases quite substantially, with masses derived from 
Eq.5 being on average 40 per cent smaller than the true ones. 
The discrepancy between true and recovered cluster masses is 
larger than in B04 as a consequence of using T^\, instead of T^^. 
A comparable bias in the mass estimates has been also found 
by Bartelmann & Steinmetz (1996), who used a spectroscopic 
definition of temperatures. Fitting the simulation results to Eq.4 
gives a = 1 .53 ±0.05 and log(Mo//i"'M0) = 13.40±0.02, while 
log(Mo//i"'M0) = 13.41 ±0.13 when fixing a = 1.5 (we note 
that the change of Mq is larger than the 40 per cent average de- 
crease of the temperature, due to the non-Gaussian distribution 
of the scatter in the left panel of Fig. 2). 

As already mentioned, the M-T relation is one of the key in- 
gredients in the recipes to extract cosmological parameters from 
the cluster X-ray luminosity (XLF) and temperature functions 
(XTF): for a fixed value of fiom, a larger Mq implies a larger 



mass for a fixed temperature and, therefore, a higher normal- 
ization of the power spectrum (e.g., Borgani et al. 2001; Seljak 
2002; Pierpaoli et al. 2003; Henry 2004). Huterer & White 
(2002) suggested an approximation for the scaling of cg with 
Mq, involving the matter density parameter: l^o^^crg oc Mq-^^ . 
Based on this relation, we expect that a bias in the mass esti- 
mate as large as that found for the simulated clusters turns into 
an underestimate of by about 20 per cent. 

Besides decreasing Mo, estimating cluster masses from Eq.5 
also significantly reduces the scatter in the M-T relation from 
about 30 per cent to about 16 per cent, as evident from Figure 
2. This witnesses that the complex structure of the ICM en- 
ters in determining the intrinsic scatter of the M-T relation in 
a more subtle way that just accounted for by cluster-to-cluster 
variations of the /3-model fitting parameters and of the effective 
polytropic index. The intrinsic scatter of the M-T relation is a 
further piece of information entering into the determination of 
(Tg: since the theoretical mass function needs to be convolved 
with this scatter, the high-temperature tail of the model XTF 
grows with the scatter, at a fixed erg (e.g. Pierpaoli et al. 2003). 
As a consequence, an overestimate of the scatter turns into an 
underestimate of ug. We find that erg decreases by about 5 per 
cent by assuming 30 per cent, instead of 16 per cent, intrinsic 
scatter Although this effect only partially compensates for the 
underestimate of the amplitude of the M-T relation, it high- 
lights once again how precision measurements of cosmological 
parameters require having under control a number of system- 
atic effects entering in the conversion between cluster mass and 
observable quantities. 

As a word of caution, we remind that the mismatch between 
Tew and T^i is introduced by the complexity of the cluster ther- 
mal structure and, therefore, so does the estimate of any biasing 
in the mass estimate. Furthermore, in order to properly quantify 
the amount of a possible bias in the observational determination 
of cluster masses it is required to reproduce as closely as pos- 
sible the observational procedure and the scale range where the 
value of /3fit is fitted (see the discussion in B04). We desert 
to a forthcoming paper (Rasia et al., in preparation) a detailed 
analysis of the observational systematics in the determination of 
cluster masses from X-ray data. In this respect, it will be also 
important to compare results on erg from the XTF with those 
from other methods, also based on galaxy clusters, but not rely- 
ing on the choice of the M-T relation. An example is provided 
by the approach based on the cluster baryon mass function (Vo- 
evodkin & Vikhlinin 2004), which provides erg ~ 0.7. 

In Figure 3 we show the effect that using T^i, instead of Tew, 
has on the estimate of the XTF. We compare here the cumula- 
tive temperature function A^(> T) from the simulation (Set#l 
only) to the data by Ikebe et al. (2002). Using Tg„ reproduces 
the observed XTF quite well (see also B04): this would sup- 
port cTg = 0.8, the value assumed in the simulation. However, 
when Tsi is used to better account for the observational proce- 
dure of temperature measurement, the simulation XTF is sig- 
nificantly affected, especially in the high-temperature tail. This 
result shows that erg = 0.8 is excluded at 3cr level, while favor- 
ing a larger power spectrum normalization, at least for the ICM 
physics introduced in our simulation. For instance, at 2 keV the 
simulated cumulative XTF is lower than the observed one by 
about 50-60 per cent. The value of Mvii corresponding to this 
temperature in our simulation is about 1.5 x lO''*/i~'M0. Com- 
puting the mass function by Sheth & Tormen (1999) for our 
cosmology at this mass, we find that erg needs to be increased 
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Fig. 3. — The cumulative temperature function A'(> T), from Set#l. Filled 
squares and dashed line refer to the results for the emission-weighted temper- 
ature Tew, while filled triangles and solid line represent the spectroscopic-like 
temperature Ts\. Error bars correspond to Poissonian uncertainties. For refer- 
ence, open circles show the local temperature function measured by Ikebe et 
al. (2002) and adapted for the considered cosmological model. 



to about 0.9 to increase the number density of clusters above 
this mass hmit by the required 50-60 per cent. 

4. CONCLUSIONS 

As demonstrated by Mazzotta et al. (2004; see also Math- 
iesen & Evrard 2001), the emission-weighted temperature, Tg„, 
usually assumed to analyze the results of hydrodynamical sim- 
ulations, fails in reproducing the "observational" temperature 
obtained from the spectral fit. As an alternative, M04 proposed 
and tested a new formula, named spectroscopic-like tempera- 
ture, Tsi, which gives a much better approximation to the spec- 
troscopic temperature, with differences of a few per cent for 
T > 2-3 keV. 

In this Letter we have studied the consequences, for the mass- 
temperature (M-T) relation and for the estimate of cts from the 
X-ray temperature function (XTF), of adopting Tsi, instead of 
Tew We have analyzed a set of hydrodynamical simulations of 
galaxy clusters, which include a realistic treatment of the gas 
physics. Our main results can be summarized as follows. 

(a) We find in our simulations that T^i is smaller than Tew by a 
factor 0.7-0.8. A linear fit approximating the relation between 
the two different temperature estimators is given by Eq. 3. 

(b) Using Tsi, instead of Tew, increases the normalization of 
the M-T relation making it higher than the observed one by 
about 50 per cent. By assuming hydrostatic equilibrium for a 
gas density distribution described by a /3-model with a poly- 
tropic equation of state, masses are underestimated on average 
by ~ 40 per cent, therefore substantially reducing the discrep- 
ancy with observational data; the scatter in the M-T relation is 
also underestimated by about a factor two. 

(c) A bias in the M-T relation propagates into a bias in the 
power spectrum normalization, a^, from the XTF. If such a bias 
is as large as that found in the simulations, the values of erg 
obtained by combining the local XTF and the observed M-T 
relation are underestimated by about 15 per cent. 

(d) The XTF from the simulation is significantly lower when 



using Tsi instead of Te„. A comparison with the observed XTF 
(Ikebe et al. 2001) indicates that for the "concordance" ACDM 
model (78 needs to be increased from ~ 0.8 to ~ 0.9. 

In general, our conclusions go in the direction of alleviat- 
ing a possible tension between the power-spectrum normahza- 
tion obtained from the number density of galaxy clusters and 
that arising from the first-year WMAP CMB anisotropics (e.g. 
Bennett et al. 2003) and SDSS galaxy power spectrum (e.g. 
Tegmark et al. 2004). Our results rely on: /) the reliability of 
our simulations to correctly describe the thermal complexity of 
the ICM; //) the capability of reproducing and understanding 
any possible source of bias in the estimates of cluster masses. 
There is little doubt that these two issues need to be addressed 
in detail in order to fully exploit the role that hydrodynamical 
simulations play in the calibration of galaxy clusters as preci- 
sion tools for cosmology. 
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